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The standard model for describing propagation of a pulse-shaped complex field 
envelope in nonlinear dispersive media is the nonlinear Schrodinger (NLS) equation. 
In the context of nonlinear optics, the main assumption made when deriving the NLS 
equation from Maxwell's equations is that the pulse-width is large as compared to the 
period of the carrier frequency. When this assumption is no longer valid, i.e., for pulse 
duration of the order of a few cycles of the carrier, the evolution of such "short pulses" 
is better described by the so-called short-pulse equation (SPE) pp. 

The SPE can be expressed in the following dimensionless form, 

u xt = u + \{v?) xx (1) 
o 

where subscripts denote partial derivatives. The SPE forms an initial-boundary problem 
when accompanied by the initial data 

u(0,x) = u Q , 

and sufficiently fast decaying boundary conditions u(t, ±00) = 0. Much like the NLS 
equation, the SPE is integrable [2] and exhibits soliton solutions in the form of loop- 
solitons [3J. However, when it is formed as an evolution equation certain conditions 
must apply otherwise, as shown below, the temporal derivative exhibits discontinuities. 

Despite the fact that the equation is integrable via the inverse scattering transform 
[I], there are certain subtleties that need to be clarified. Integration of Eq. ([I]) introduces 
the operation 

/X 
u(t, x') dx' 
■00 

Clearly as x approaches —00, d~ l u = 0, consistent with rapidly decaying data. However, 
as x approaches +00, for u and its time and space derivatives to decay, a constraint 
seems to be necessary (see the discussion below), namely 

X) 

u(t, x)dx = (2) 

> 

Indeed, writing the SPE in evolution type form we have 

1 f x 1 

u t = d~ l u + = / udx'+-(u% 

6 J-00 6 

and imposing the boundary condition as x — > +00, one results to Eq. (j2J). 

In fact, this constraint induces further constraints obtained by successively taking 
the time derivative of the integral and using Eq. ([T|). For example, the next constraint 
is given by 

p+00 

d~ x u dx = (3) 

> 

However, Eqs. (J2J), (J3J, along with the rest of the family of infinite constraints generated 
as above, are not generically true. One might surmise that constraints are required at 
all times for a solution to exist. However, as discussed below, this is not the case. 
Extra constraints on the initial data are not necessary, but the solution suffers from 
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a temporal discontinuity. For smooth initial data not satisfying Eq. (j2j), u t (t,x) has 
at t — different left and right limits and the rest of the family of constraints cannot 
be generated dynamically at that point. The same issues arise in the context of the 
Kadomtsev-Petviashvili (KP) equations and were studied in Refs. [SI [6]. 
Our analysis starts by taking the Fourier transform (FT) of Eq. (CQ), 

k 2 - 

iku t = u — — u 3 (4) 
6 

where the FT pair is defined as 

/+oo 
u(t,x) e ikx dx 
■oo 
i r+oo 

u(t, x) = T~ l {u(t, k)} =— u(t, k) e~ ikx dk 
Define U = u 3 and write Eq. (jl]) in the form of a first order differential equation in t, 

which can be readily integrated with the use of integrating factors to give 

u(t, k) = u e t/ik + — e t/ik f U e~ T/ik dr (6) 



o 



6 

where u = Uo(k) = JF{m (x)}, is the FT of the initial data. Using Eq. we calculate 
the temporal derivative to be 



(t, k) = |n e t/lfc + ~7rU + ~ flJe^ k dr (7) 
ik 6 6 J 

Clearly, from Eq. (jSJ), as tends to zero, we should demand that so will u. This 
translates to [7] 

/■+O0 

u(t,0) = <=> / u(t,x) dx = 



However, as t — > ±0, we have that u(0, x) = Uq from Eq. ([6]), and 

1 . ik* 

Ut = 1 : — ?7vT u ° + T U 

ik — sign(rJL) o 

This is because the function exp(t/ik) defines a distribution, depending continuously on 
t, in the Schwartz space of the variable k jSj with 



and 



^ ttlk = - 1 -——e u / k , t = 

at ik — sign(t)0 



at ik 



This suggests that although there is no discontinuity in the solution, there is one in the 
derivative. Indeed, taking the inverse FT of Eq. ([7]), at t — >• ±0, we have 

i r+°° i i 

u t (t -> ±0,x) = — lim / — u e t/ik e~ ikx dk + -(u 3 ) x (t -> ±0,x) 
2ti t->±o J_ 00 ik 6 
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The nonlinear term is straightforward to handle so we focus on the linear part, 



Kx) 



— lim 

27T t-*±0 



oo j 

— u e t/ik e~ ikx dk 

IK 



oo 
' oo 



1 f + °° 1 

2^ So ik ^~ lkX + Uo{0) ~ ^ dk 



— lim 

27T t~>±0 



— lim 

27T t— ±0 



OO 

oo 



CX) 



/A' 



[M e" ifc3; - tio(O)] e t/ik dk 



^ru (0)e t/ik dk 
ik 



(8) 



Using the property 



— e t/ik dk = -7r sign(t) 

„ IK 



the second integral of Eq. (jSJ) is reduced to — uo{0)tt sign(£)/2. Furthermore, we write 



u (0) 



6(k)u (k) e~ lkx dk 



so that finally 



Kx) 



l 

2^ 
1 

2^ 



— oo 
oo 



1 

ik 
u (k) 



ik — sign(t) 
Uo(x') dx' 



7r sign(t)5(/c) 



M (fc) e" 



dA; 



'sign(t)oo 

where P denotes principal value. Thus, at i = 0, Eq. <^7§ translates into physical space 
as 



Ut 



sign(t)oo 



Uq(x') dx' + ~(Uq) 
o 



As also mentioned in Refs. [HE], the operator 



d- 1 



dx' 



sign(t)oo 



and its relative average 



dx' + 



dx' 



are equivalent, meaning that one can choose either one of them. If t ^ we have that 
u dx = 0, hence both choices are valid. At t = there is a discontinuity in the 
temporal derivative. 

For the evolution of the SPE, Eq. ([2]) is not preserved in time and as such leads 
to the infinite number of further constraints. Indeed, if Eq. (TS]) holds then an infinite 
number of constraints, dynamically generated using the SPE, hold during the evolution. 
Within the physical framework of the SPE these constraints are neither "natural" nor 
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necessary. Solutions of the SPE can exist without satisfying this condition, the most 
prominent example being the loop-soliton [3] . This solution, however, in addition to the 
possible temporal discontinuities, suffers from discontinuities in its spatial derivatives, 
u x (t,x), and extra care may be needed when the above formalism is applied. 

We conclude with a note on the so-called regularized SPE (RSPE) model, recently 
derived in Ref. [9j. The latter has been derived by including a regularization term, 
based on the next term in the expansion of the dielectric's susceptibility. In that case, 
the pulses (of the real component of the electric field) are described by: 

O 

where (3 is a small parameter. Without the regularization term, /3u xxxx , i.e., in the case 
of the SPE -cf. Eq. ([I])-, traveling pulses in the class of piecewise smooth functions with 
one discontinuity do not exist. However, when the regularization term is added, and for 
a particular parameter regime, the RSPE supports smooth traveling waves which have 
structure similar to solitary waves of the modified KdV equation [9]. The regularization 
term does not alter the analysis for the SPE. Indeed, in the Fourier domain the term is 
written as T {f3u xxxx } = f3k A u and thus when dividing with ik from the left-hand-side 
the resulting power of k is continuous at k — 0. The linear part of the RSPE, much like 
the linear part of the KP-I equation [51 [8] , deserves more study and the analysis will be 
presented in a future communication. 
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